Energy-sensitive imaging detector applied to the dissociative recombination of D 2 H 
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We report on an energy-sensitive imaging detector for studying the fragmentation of polyatomic 
molecules in the dissociative recombination of fast molecular ions with electrons. The system is based 
on a large area (10x10 cm 2 ) position-sensitive, double-sided Si-strip detector with 128 horizontal 
and 128 vertical strips, whose pulse height information is read out individually. The setup allows 
to uniquely identify fragment masses and is thus capable of measuring branching ratios between 
different fragmentation channels, kinetic energy releases, as well as breakup geometries, as a function 
of the relative ion-electron energy. The properties of the detection system, which has been installed 
at the TSR storage ring facility of the Max-Planck Institute for Nuclear Physics in Heidelberg, is 
illustrated by an investigation of the dissociative recombination of the deuterated triatomic hydrogen 
cation D2H + . A huge isotope effect is observed when comparing the relative branching ratio between 
the D 2 +H and the HD+D channel; the ratio 2B(D 2 ±H)/B(HD±D), which is measured to be 
1.27 ±0.05 at relative electron-ion energies around eV, is found to increase to 3.7 ±0.5 at ~ 5 eV. 

PACS numbers: 34.80. Lx, 34.80.Ht, 34.80.Gs 



I. INTRODUCTION 

Dissociative recombination (DR) with electrons is the 
most important loss channel for molecular ions in cold di- 
lute plasmas and a source of energetic neutral atoms and 
excited molecules pj. Rate coefficients, branching ratios, 
and the energy sharing between internal and external de- 
grees of freedom of the fragments are thus important in- 
gredients in modeling these plasmas. Even though a good 
understanding of the DR process of diatomic molecular 
ions has been reached in recent years, theory is not yet 
capable of supplying most of the information needed for 
these studies. Experimental investigations of the DR pro- 
cess in particular of polyatomic systems are thus still in- 
dispensable to further constrain the parameters of these 
models and to provide basic theory with detailed and 
reliable results for crucial benchmark systems. 

Considerable experimental progress in studying the 
DR process has been made since the advent of heavy-ion 
storage rings with merging electron beam facilities, which 
allow to control the inner excitation of the molecular ions 
and the relative ion-electron energy precisely enough to 
produce data in the parameter space of interest for, e.g., 
planetary atmospheres or interstellar chemistry. Results 
obtained so far were vital in helping to unravel the pro- 
cesses governing the DR and to clarify many important 
and pertinent issues. But despite these experimental ad- 
vances for particular diatomic systems, detailed investi- 
gations of the DR of more complex ions, such as measure- 
ments of branching ratios at relative energies other than 
zero, of the inner excitation of the molecular fragments or 
of their dissociation kinematics, are still hampered due to 
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the limitations of presently available detection techniques 
in identifying the DR reaction products. 

We have therefore developed an Energy-sensitive 
MUlti-strip detector system (EMU) that is capable of 
recording multi-fragment events following the DR of 
polyatomic molecular ions and to identify the fragments 
by their masses. The system has been installed at the 
ion storage ring TSR of the Max-Planck Institute for 
Nuclear Physics, which can store molecular ions of McV 
energies and is equipped with an electron cooler as well 
as an electron target employing a cryogenic photocath- 
ode, so that electron-ion collision energies can be freely 
set in the mcV to cV range with unprecedented resolu- 
tion Q. While the main aim of the setup is to provide 
a universal tool for measuring DR branching ratios of 
polyatomic molecular ions into the different final frag- 
ment channels, the position resolution is sufficient to al- 
low also for 2D imaging of the breakup geometries. The 
new setup has thus several advantages as compared to 
previous approaches Q , where grids of different transmis- 
sions (yielding a set of linear combinations of the rates 
for individual fragment channels) were employed in or- 
der to determine the branching ratios, a technique which 
is usually only applied at zero relative electron-ion ener- 
gies for background reasons and which does not support 
molecular imaging. 

The present paper discusses the concept and realiza- 
tion of the EMU detection system and exemplifies its 
properties and possibilities by using the new setup to in- 
vestigate the dissociative recombination of D2H + . The 
DR of D2H + was studied before in several storage ring ex- 
periments, and rate coefficients [1, Q , the fragmentation 
geometry of the three-body channel @, and branching 
ratios [5| have been reprted. The kinetic energy Qo re- 
leased in the three fragmentation channels, for ions and 
fragments in their ground states and relative ion-electron 
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energies E e = eV, arc 

( H + D + D (+4.67 eV) (a) 
D 2 H + +e" — >\ HD + D (+9.19eV) (/3) . (1) 
[ D 2 +H (+9.23 eV) ( 7 ) 

The branching ratios measured at E e ~ eV were 
reported to be 5 Q =76.5(2.2) %, 5^=13.5(1.5)%, and 
7^=10.0(0.7)% @. 

II. THE DETECTOR CONCEPT 

Since the beam energies E B of the molecular ions 
stored in the TSR are very large as compared to the 
release energies occurring in DR experiments, the frag- 
ments resulting from a DR event of a molecular ion 
of mass number A are traveling within a narrow cone 
around the ion beam axis, and the kinetic energy E^ of 
the fragment i is to a very good approximation given by 

Ei = ^jE B (2) 

and is thus a measure of the mass number Aj of the 
fragment. By recording separately the energies of all 
DR fragments of an DR event with a large area multi- 
hit-capable detector the fragmentation channel can be 
uniquely identified. 

In the present setup the energy-sensitive multi-hit de- 
tector is realized by a large area (10x10 cm 2 ) position- 
sensitive, double-sided Si-strip detector (see Fig.Q}. The 
energy- and position-sensitivity to particles hitting the 
detector is achieved through 128 vertical (x-) and 128 
horizontal (y-) strips on the front and back side of the 
detector, respectively, which are read out individually by 
preserving the pulse height, i.e., the energy information. 
The readout system also ensures via the timing informa- 
tion that the detected particles belong to a single DR 
event. The detector is mounted about S sw 10 m down- 
stream of the electron target outside of the magnetic lat- 
tice of the TSR such that it can only be hit by neutral 
fragments. 

The impact of n fragments from a single DR event on 
the detector leads to n x responding x-strips Xi and n y 
responding y-strips yi with corresponding energies E Xi 
and E Vi , and the detection of all fragments can be easily 
verified by checking whether the two energy relations 

g E Xi = E B and E v< = E b (3) 

i=l i=l 

are fulfilled. The measured fragment mass numbers A Xi 
and A yi corresponding to the energies recorded on the 
individual strips are given by 

A Xi = AE Xi /E B and A Vi = AE yi /E B . (4) 

In the following we will refer to the set of measured mass 
numbers A Xi , A yi and the corresponding strip coordi- 
nates Xi, yi as the 'hit pattern' of the DR event. 



The coordinates of the center-of-mass (x cm , y cm ) of the 
detected event in the detector plane are given by 

/ Tlx n v \ 

{XcrmVcm) = l/-^ I ^ ^^ A Xi Xj, ^ A yi yj J , (5) 
\i=l i=l / 

and an estimate for the transversal kinetic energy release 
E± can be determined from the projected weighted dis- 
tances D 2 defined by 

Tlx a n y * 

V 2 = E ^r(*< - *-) 2 + E ^(** - Vcrnf (6) 
i=l i=l 

using 

E ± =E B x D 2 /S 2 , (7) 

where S is the distance between the center of the electron 
target and the detector. 

While Eqs. ©-([7]) are valid even in cases where more 
than one fragment hit the same strip, the strip- wise read- 
out of only two coordinate planes together with the fi- 
nite strip width leads to ambiguities in the determination 
of the breakup geometry from the hit pattern, which in 
some cases also affects the identification of the fragmen- 
tation channel. We shall first discuss these ambiguities 
for the specific case of the DR of D2H" 1 ", before we con- 
sider the more general case. 

Figure [2] depicts the two most frequent of these ambi- 
guities. The hit pattern shown in Fig. [^a) is indistin- 
guishable from the hit pattern shown in Fig. [2jb); this 
pattern thus results in two solutions concerning the frag- 
mentation geometry. However, the identification of the 
fragmentation channel is still unique. This is no longer 
the case for the hit pattern shown in Fig.[2jc), where the 
H and the two Ds from the 3-body fragmentation channel 
are forming an L-shaped pattern with the H in the corner 
such that only two horizontal and two vertical strips are 
responding. This pattern cannot be distinguished from 




S~10 m 



FIG. 1: Concept of the EMU detector: A fast, neutralized 
molecule is dissociating at a distance S from a position- and 
energy-sensitive detector. The fragments continue to travel 
approximately at the velocity of the molecule, but due to 
the kinetic energy released in the dissociation they acquire 
macroscopic distances from the transversal center-of-mass 
(x C m,ycm) when hitting the detector. The position (xi,yi) 
and the energy Ei of each fragment is recorded. 
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FIG. 3: Schematic cross section through the EMU detector, 
a double-sided, fully depleted Si-strip detector of 10 x 10 cm 2 
active surface area. The p + Si and the Al layers of the junction 
contact, through which the particles enter the detector, were 
made as thin as possible to minimize the energy losses of the 
DR fragments in these dead layers. 



FIG. 2: Some breakup geometries and corresponding hit pat- 
tern occurring in the DR of D2H + . The lines mark the re- 
sponding read-out strips in vertical (xi) and horizontal (yi) 
directions; the numbers at the end are the mass numbers A Xi 
and A Vi derived with Eq. (j4| . The breakup geometries shown 
in a) and c) result in the same read-out pattern (hit pat- 
tern) as the geometries shown in b) and d), respectively. But 
while the hit pattern resulting from the geometries shown in 
a) and b) still allow unique identification of the fragmenta- 
tion channel, the pattern resulting from the breakup geome- 
tries depicted in c) and d) is ambiguous with respect to the 
fragmentation channels D+D+H and HD+D. 



(sub-)pattern involving two identical fragments at the 
two ends of the L. The interpretation of hit patterns in 
terms of the fragmentation geometries is subject to sim- 
ilar ambiguities as discussed above, but the number of 
geometries leading to a specific pattern will strongly in- 
crease with the number of identical fragments. However, 
we again would like to point out that the center-of-mass 
determination by Eq. ([5]) and the transverse energy dis- 
tribution defined by Eq. (J7J arc independent of these am- 
biguities. 



that shown in Fig.^d), which results from the detection 
of a molecular HD- and an atomic D-fragment from the 2- 
body channel (/3). Besides the even less likely occurrence 
of two or more fragments hitting the same x- and y-strip, 
i.e., the same (x, y)-pixel, these L-shaped patterns are 
the only patterns in the DR of D2H + , which cannot be 
attributed uniquely to a fragmentation channel. Since 
the probability for the occurrence of this pattern is small 
due to the narrow strip width of 730 /nm realized in the 
present setup, we attribute these patterns to the channel 
with the smaller number of fragments, that is for the case 
discussed in Figs. (2{c,d) to the HD+D channel. The in- 
fluence of these miss-assignments on the accuracy of the 
branching ratios is small and can moreover be corrected 
for with the help of Monte Carlo simulations (see also 
Sec. lVA]) . 

While the ambiguities caused by fragments hitting the 
same (x, y) pixel discussed above are present in case of 
all molecules, the other ambiguities depend on the mul- 
tiplicity of identical fragment masses in the open frag- 
mentation channels. In the DR of polyatomic molecules 
consisting of atoms of different masses such as HCO + , 
the identification of the fragmentation channel and of 
the fragmentation geometry is not subject to any other 
ambiguities. In the DR of polyatomic molecules, which 
contain several atoms of the same mass and fragment 
into channels containing two or more identical fragments, 
the main ambiguities hampering the identification of the 
fragmentation channels arc again due to the L-shapcd 



III. EXPERIMENTAL SETUP 
A. The EMU detection system 

The large area double-sided Si strip detector employed 
in the EMU detection system was built by the UK physics 
company Micron Semiconductor Limited 0]. The detec- 
tor has an active area of 97.3 mmx 97.3mm and an ac- 
tive depth of 300 /im when biased by < —70 V (Fig. [3]). 
Position sensitivity has been achieved by electronically 
subdividing the junction side into 128 vertical strips of 
730 /im width, which are separated by gaps of 30 /im, and 
the ohmic side in 128 horizontal strips of 700 /^m width, 
which arc separated by gaps of 60 /im. 

The thickness of the entrance (junction) window, 
which consists of a lOOnm-thick p + doped Si layer cov- 
ered with a lOOnm-thick aluminum coating, was mini- 
mized in order to keep the energy loss in the inert window 
material, and thus the cut-off energy for the detection of 
the heavy DR fragments, as small as possible. The in- 
creased thicknesses of the p + Si and the Al layer at the 
edges of the strips, which ensure a clean separation be- 
tween the strips and a loss-free readout and which affect 
about 14 % of the active area of the detector, result in ad- 
ditional, slightly down-shifted peaks in the pulse height 
spectra; particles traversing this part of the entrance win- 
dow suffer a larger energy loss which is specific for the 
fragment energy, mass, and its nuclear charge. This is 
discussed in more detail in Sec. IIVI in connection with 
the pulse height spectrum (Fig. [5]) observed in the DR of 
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FIG. 4: Schematic view of the EMU setup, which is located 
in the BAMBI beamline downstream of the electron target of 
the TSR. 



to do all the controlling, timing and read out of the 16 
MTM-16 boards. A trigger signal is produced when one 
of the discriminators of the 16 MTM-16 units responded; 
it activates the hold stages of the MTM-16 units after 
2 fxs, (i.e., only signals arriving within this time will be 
held) and starts the readout sequence. The individual 
pulse heights are digitized by a sliding scale ADC and 
- if above a programmable threshold - stored together 
with the strip number in a FIFO, which is read out via 
the VME interface. In the present configuration the read- 
out speed is limited to about 2000 events per second. A 
schematic view of the electronic set up is shown in Fig. [5] 



D2H + . Since the penetration depth of the fragments in 
the detector material for energies available at the TSR is 
in most cases less than 10 fim and the gap region close 
to the surface is almost field-free, fragments impinging 
in the gap region, covering about 4 % of the detector 
surface, will likely not be detected. 

The detector is mounted in a dedicated chamber 
(Fig. [3J in the neutral fragment beamline BAMBI down- 
stream of the electron target of the TSR. The distance 
from the center of the electron target to the front face of 
the detector is So = 941(1) cm. At this distance the de- 
tector covers the full cone size for the DR fragments that 
is allowed by the bcamtubc (CF100). The maximum vis- 
ible cone is in fact limited vertically by the height of the 
vacuum chamber in the TSR dipolc magnet to about 8 cm 
diameter. A large bellow allows to retract the detector 
out of the beamline so that experiments with detectors 
further down the beamline are possible without break- 
ing the vacuum. Moreover, the detector can be moved 
up and down in front of a retractable a-source, which 
in turn can be moved horizontally using a translational 
stage, so that any part of the detector can be irradi- 
ated for testing and calibration purposes. The cryopump 
is vibrationally decoupled from the detector chamber to 
protect the wire bonds to the 256 strips of the detector 
from the mechanical noise of the pump. Under running 
conditions the vacuum in the detection chamber is about 
5 x 10 -9 mbar. A mechanical shutter is furthermore in- 
stalled in the pumping chamber to protect the detector 
during the injection phase of the ion beam into the stor- 
age ring. 

The 256 strips of the detector are read out by 16 highly 
integrated MTM-16 units developed by Mesytec y|, 
which arc mounted directly on the upper lid of the detec- 
tor chamber to minimize the length of the cables between 
detector and preamplifier. Each unit contains 16 chan- 
nels consisting of a preamplifier, main amplifier, shaper 
and hold stage. Moreover, four input signals are added 
and compared to an adjustable discriminator level. The 
eight boards connected to the 128 x-strips and the eight 
boards serving the y-strips are daisy-chained to the two 
function blocks, respectively, of an MDI-2 unit developed 
by Mesytec as well. This one-slot VME unit allows one 



B. Specific settings for measuring the DR of D2H + 

The first application of the EMU detection system 
was the investigation of the dissociative recombination 
of D2H+. D2H + ions were produced in a Penning source, 
accelerated to Eb = 2.39 MeV using an rf-quadrupole ac- 
celerator, and injected into the TSR; the stored ion cur- 
rent was up to about 4 fj,A. After injection the ion beam 
was first phase-space cooled by the velocity-matched elec- 
tron beam supplied by the electron target, which had a 
density of 9xl0 6 cm~ 3 and a transversal temperature of 
approximately 1.5 meV (expansion factor 20). At these 
conditions the effective length of the electron target is 
115 cm, with an additional 20 cm on either end where the 
electron beam is merged with and demerged from the ion 
beam. After a few seconds the ion beam was phase-space 
cooled to a beam diameter below 1.0 mm. 

The electron beam was then used as a target, cither 
at matched or detuned electron velocities. The detuning 
energy Ed, which measures the average electron-ion col- 
lision energy in the center-of-mass system, was adjusted 
to values between and 20 eV by changing the accelera- 
tion voltage of the electron beam. At detuning energies 
Ed > the electron beam velocity was switched at a rate 
of 10 - 100 Hz between Ed = and the required energy 
Ed in order to keep the ion beam phase-spaced cooled. 
Data was recorded between a few seconds and up to 30 s 
after injection, and count rates were limited to < 2 kHz. 
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FIG. 5: A schematic view of the electronic setup for the EMU 
detector. 
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IV. DATA ANALYSIS 

The first analysis step consists of shifting the pulse 
height spectra of the individual strips to a common scale 
in order to compensate for the slightly different amplifi- 
cations and offsets of the individual electronic channels 
using the D 2 H+ data. 

The sums of all 128 x- and of all 128 y-strip spec- 
tra performed after this calibration step, are displayed 
in Fig. |51 While the overall energy resolution reached 
with the y-strips on the back side of the detector (e.g., 
A£(FWHM)< 30 kcV for Ds of 0.96 MeV) is already 
very satisfactory, the resolution of the front side strips is 
worse by almost a factor of two. This is mainly due to the 
excess noise observed on the front side strips; while the 
back side of the detector is electronically well-shielded 
by a metal plate mounted directly behind the detector, 
the front side is unshielded in the direction of the storage 
ring. But although the resolution of the x-strips is not 
yet perfect, the peaks corresponding to the detection of 
an H, D, HD (or H+D), and D 2 (or D+D) can be clearly 
identified. Even the small peak is visible that is caused 
by the detection of all fragments (D2H) on a single strip. 

On the low-energy side of the peaks small satellite 
lines are visible. As already pointed out in Sec. IIII Al 
these satellites are caused by particles passing through 
the slightly thicker parts of the entrance window. For 
light fragments as in the present case of the DR of D2H + 
the resulting energy difference between the main and the 
satellite peak is small and thus does not hamper the mass 
assignment. For heavier fragments like carbon or oxygen, 
the energy difference will be larger and can correspond 
to one or more mass units. But even in these cases this 
will not lead to a misinterpretation of a DR event as the 
two energy sums [see Eq. will not be fulfilled. Dis- 
regarding these events, however, reduces the detection 
efficiency of the set-up and has to be taken into account 
by detailed simulations. 

A more careful examination of Fig. |U[b) reveals a small 
number of counts between the peaks. These are mainly 
caused by events where the charges created by a frag- 
ment in the active volume of the detector are drifting to 
two adjacent y-strips such that the charge collected by 
each strip results in a reduced signal height. Such a pulse 
height splitting can be at least partly reconstructed by 
adding the signals of the two adjacent y-strips. A corre- 
sponding pulse height splitting on the front side of the 
detector is very unlikely to occur [see also Fig. [SJa)] as 
the range of the fragments in the detector material is so 
small that fragments hitting the gap between the x-strips 
are likely not be detected at all. As in the case discussed 
above, the influence of these effects on the detection effi- 
ciency is taken into account by simulations. 

After the calibration procedure the pulse height of each 
of the 256 strips is compared to appropriately chosen 
windows around the mass peaks and, if it falls into one 
of these windows, the corresponding mass number is as- 




200 400 600 800 1000 1200 1400 
channel number 

FIG. 6: Pulse height spectra obtained in the DR of 2.39 MeV 
D2H + ions at a detuning energy of Ed = eV. The figure 
shows the sum of the calibrated pulse height spectra (a) of all 
x-strips on the front and (b) of all y-strips on the back side. 
Above channel 550, counts are multiplied by a factor of 10. 



signed to the strip. If the pulse height is above the lower 
threshold but does not fall into one of the mass windows, 
it is tested whether the sum with one of the neighbor- 
ing strips fulfills this criterion. If yes, the corresponding 
mass is assigned to the strip with the larger signal while 
the other strip is ignored in the further analysis. We 
shall refer to this set of A Xi and A Vj and adjusted posi- 
tions Xi and yj, which are randomized across the width 
of the respective strip to avoid binning problems, as a hit 
pattern. 

The next analysis step consists of checking whether 
the sums over all mass numbers A Xi and A yj , respec- 
tively, observed in an event are equal to the mass num- 
ber of the dissociating molecule, and whether the hit pat- 
tern corresponds to a possible DR fragmentation channel, 
i.e., whether there is any combination of the atomic con- 
stituents of the dissociating molecule that can explain the 
observed hit pattern. If yes, the event is marked as a DR 
event and the final fragment channel is assigned. Here, 
preference is given to the fragmentation channel with the 
smallest number of fragments. In the case of the DR of 
D2H + , for example, this means that the rather unlikely 
L-shaped hit pattern shown in Fig. [He), which results 
from the three-body breakup and cannot be distinguished 
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FIG. 7: The distribution of the H (black) and D 2 (red) po- 
sitions observed for 5000 D2+H events following the DR of 
D2H + at Ed=0eV. Note that H events within the red area 
are partly shadowed by the D2 events due to the plotting 
procedure. 



from the hit pattern caused by the two-fragment chan- 
nel HD+D as the comparison with Fig. ^d) shows, is 
assigned to the HD+D channel. 

An example for the extracted position distribution of 
fragments for a single fragment channel is given in Fig. [3 
where the positions of the D2 (red) and the H (black) 
fragment of the D2+H channel following the DR of D2H+ 
at Ed = eV are shown in form of a scatter plot. One 
can clearly see the two circular shapes corresponding to 
the fragments with mass numbers 1 and 4, respectively. 

From the positions Xi and yj and the mass num- 
bers A Xi and A yj the (transversal) center-of-mass po- 
sition (x cm , y C m) is calculated. The resulting x- and 
y-distributions of the center-of-mass position for events 
following the DR of D2H + at Ed = OeV are shown in 
Fig. [8] The data was taken 30 s after injection and con- 
tinuous electron cooling, and exemplifies the good phase- 
space cooling that can be reached with the electron tar- 
get of the TSR; taking the width of the strips into ac- 
count the center-of-mass spread at the position of the 
detector is estimated to be < 0.9 mm (FWHM) in x- as 
well as in y-direction. In the further analysis of the data 
only (-vents within a. 2rj-cllipse around the centroid of the 
center-of-mass distribution are used to further suppress 
background events which might result from dissociative 
charge exchange reactions of the molecular ion with the 
residual gas of the TSR. 

To investigate the influence of the detector properties 
and of the approximations and cuts used in the data anal- 
ysis on the deduced results, and to allow for a detailed 
comparison with theoretical predictions for the DR pro- 




x-position (strip number) y-position (strip number) 

FIG. 8: The horizontal and vertical distribution of the center- 
of-mass coordinates of DR events of D2H + taken after 30 s of 
phase-space cooling at collision energies of Ed — OeV. The 
blue lines are fits by a Gaussian distribution. 

cess, a Monte Carlo simulation program has been devel- 
oped. The simulation creates DR events at random posi- 
tion in the overlap region between the molecular ion beam 
and the electron beam and propagates them towards the 
detector. Subroutines are available to allow for an in- 
ternal excitation of the ion and to incorporate different 
kinetic energy releases and angular distributions of the 
fragments with respect to the center-of-mass direction. 
At the detector, the impact position for each fragment is 
determined, which includes the decision as to which part 
of the entrance window was hit. The kinetic energy loss 
in the appropriate dead layer is subtracted for each frag- 
ment and a pulse proportional to the remaining kinetic 
energy is created on the strip being hit, broadening the 
pulse height by a Gaussian shaped resolution function. 
For fragments hitting the same strip the pulse heights 
are added. If the gap between two x-strips is hit, the 
fragment is assumed not to be detected. If the impact 
position happens to be in the gap between two y-strips 
of the back side, the signal is shared. The output of the 
program for a simulated event is a list of strip numbers 
and corresponding pulse heights. The simulated data are 
then treated with the same analysis tools as the measured 
data. 



V. RESULTS FOR THE DR OF D 2 H+ 

A. Branching ratio at Ed = eV 

The branching ratios for the three final fragment chan- 
nels accessible in the DR of D2H + at relative ion-electron 
energies of Ed = eV were determined by analyzing 
data taken after 30-50 s of electron cooling. Although 
the phase space cooling of the ion beam and the vibra- 
tional cooling of the ions is achieved already after ~ 3 s 
of electron cooling, the cooling was continued in order 
to lower the rotational temperature of the ions to val- 
ues below the ambient room temperature of 300 K as 
shown in Ref. The relative number of DR events 
leading to the three fragmentation channels (see Eq. ([T])) 
arc iV Q =76.69(4)%, ^=14.32(4)%, and iV 7 =8.99(3) %, 
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TABLE I: Branching ratios obtained in the present exper- 
iment for the three fragment channels following the DR of 
D 2 H+ at E d = OeV. The results are compared to those of a 
recent measurement at CRYRING [jj. 



Channel 


This work 


CRYRING [5] 


(a) H+D+D 


78.0(0.4) 


76.5(2.2) 


HD+D 


13.5(0.3) 


13.5(1.5) 


(7) D 2 +H 


8.6(0.2) 


10.0(0.7) 



where the uncertainties given are of purely statistical ori- 
gin. The Monte Carlo simulation was used to correct the 
measured numbers for the slightly different efficiencies 
for detecting and identifying the fragment channels. In 
this simulation the angular distributions of the fragments 
were assumed to be isotropic, but the population of the 
different vibrational states of the HD and D2 fragments 
as determined in Sec. IV CI was explicitly taken into ac- 
count. The resulting corrections which had to be applied 
to the measured relative numbers given above amount to 
less than two percentage points. 

The deduced branching ratios B are compiled in Ta- 
ble HI The uncertainties given are dominated by system- 
atic errors caused mainly by the not yet fully understood 
detection efficiency for fragments hitting the gap regions. 
The ratios are compared to the result of a recent measure- 
ment performed at CRYRING @ using the grid method. 
Within the error bars the two measurements agree well 
despite the presumably different initial rotational tem- 
peratures of the D2H+ ions (see also the following sec- 
tion). As discussed already in Ref. [f| the branching ra- 
tios display a clear isotope effect: On a purely statistical 
ground the HD+D channel should be twice as strong as 
the D2+H channel, i.e., in the absence of an isotope effect 
one would expect 25(D2+H)/B(HD+D) = 1. However, 
the measured values result in a deuteration enhancement 
of 2B(D 2 +H)/B(HD+D)= 1.27(0.05), slightly smaller 
than deduced in Ref. |5| but very similar to the ratio 
B(HD+HV/2B(H 2 + D) = 1.20(0.05) observed in the DR 
of H 2 D + [3(. These results clearly show that in the DR of 
deuterated H3" the formation of the two-body fragment 
channels containing the most deuterons is enhanced. 



B. Branching ratios at Ed > eV 

To determine the branching ratios into the three frag- 
ment channels at relative ion-electron energies Ed > 
eV, data collection was started after 3 seconds of elec- 
tron cooling when the ion beam was already fully phase 
space cooled. Moreover, at this time all ions are vibra- 
tionally relaxed while the rotational temperature is ex- 
pected to be around 300 K or even slightly higher jj. 
Comparing the branching ratios at Ed = OeV obtained 
after the shorter phase space cooling with those discussed 
in Sec. IV Al no difference was noticable within the lim- 
its of precision. After electron cooling, the energy of the 
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FIG. 9: (a) Energy dependence of the DR rate coefficient 
of D2H + as measured by Lammich et al. [j. (b) Fragment 
channel branching ratios determined in the present work as 
a function of the collision energy with only statistical uncer- 
tainties. In the absence of isotope effects one expects to find 
B(HD+D) = 2B(D 2 +H). 



electron beam was frequently switched between the cool- 
ing energy and the desired detuning energy Ed > eV, 
changing the latter from injection to injection to span 
the range between up to 20 eV. The detuning-energy 
cycle was repeated several times to collect statistically 
reliable data. The branching ratio at each detuning en- 
ergy was then obtained as described in Sec. IV Al includ- 
ing the corrections due to the slightly different detec- 
tion efficiencies for the three fragment channels. Effi- 
ciency losses that could occur at high detuning energies 
due to the finite size of the detector were shown to be 
insignificant by investigating the distribution of impact 
positions of the fragments on the detector; at high de- 
tuning energies a considerable part of the maximal avail- 
able release energy is not converted into kinetic energy 
of the fragments, but into internal excitation energy of 
the fragments. Moreover, background events from disso- 
ciative charge exchange collisions with the residual gas 
were found to be negligible at all measured ion-electron 
energies. 

Figure [5] shows the extracted fragment channel branch- 
ing ratios in the range of Ed = 0.2 meV to 20 eV together 
with the DR rate coefficient measured in a previous ex- 
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pcriment at the TSR under similar cooling conditions |Jj . 
As to be expected, the probability for the recombination 
channel D2H, which requires the emission of a stabilizing 
photon in times comparable to the dissociation times, is 
found to be consistent with zero to within better than 
1%. The branching ratios into the fragmentation chan- 
nels are observed to be more or less constant up to Ed « 
300 meV at values close to those obtained at Ed = OeV, 
even though the DR rate coefficient drops by two or- 
ders of magnitude. In particular, the isotope effect ob- 
served at Ed = eV persists over this energy range and 
is even increasing for energies > 300 meV. For higher 
energies the branching into the three-body channel is de- 
clining quickly and reaches a minimum level of about 
20 % between 4 and 5 eV, while the relative intensities of 
the two-body channels are both rising. At £^=2 eV the 
HD+D channel is the dominant fragmentation channel 
(B(HD+D)-45 %), but at around 5 eV the D 2 +H chan- 
nel is even reaching _B(D2+H) ~ 55 %, exhibiting a huge 
isotope effect of 2B(D 2 +H)/B(HD+D)= 3.7(0.5). At 
energies above 10 eV the three-body channel regains its 
role as the dominant fragment channel with a branching 
ratio above 60%. 

We are not aware of earlier measurements of the frag- 
ment branching ratios following the DR of D2H 4 " at de- 
tuning energies Ed > eV. However, Datz et al. 0, 
have measured these ratios in the energy range between 
~1 meV and 20 eV for the DR of H^~ and its isotopomcr 
H2D" 1 ". The overall behavior of the relative branching 
ratio between the combined two-body channel and the 
three-body channel is rather similar for all three sys- 
tems investigated, and the gross structure of this ratio 
has been successfully attributed to the successive opening 
of the electronically excited two-body fragment channel 
H(nZ)+H 2 (X%), H(20+H 2 (B 1 S„) and of the three- 
body channel H(2Z)+H(ls)+H(ls) within a statistical 
approach However, while in the DR of H2D+ the 

relative fractions of the two two-body channels seem to 
stay constant over the measured range of energies, we find 
that this fraction strongly changes in the DR of D 2 H + at 
energies between ~ 1 — 10 eV, leading to a drastic isotope 
effect around Ed ~ 5 eV that is yet to be explained. 



C. Vibrational excitation of fragment molecules 

Standard imaging techniques without mass identifi- 
cation have difficulties identifying and separating two- 
body fragmentation channels in DR experiments involv- 
ing polyatomic molecular ions. The event- by- event mass 
identification of all fragments provided by the EMU sys- 
tem allows for a straight-forward analysis of these chan- 
nels. 

The normalized projected distance distributions P(d) 
for the two two-body fragment channels D2+H and 
HD+D observed in the DR of D 2 H+ at E d = eV are 
displayed in Fig. [TU] using data collected after 30 s of 
electron cooling. Following common practice (see e.g., 




Projected Fragment Distance d (mm) 




Projected Fragment Distance d (mm) 

FIG. 10: Normalized projected distance distributions for 
the two-fragment channels D2+H (upper panel) and HD+D 
(lower panel) observed in the DR of D2H + at Ed — eV; 
the solid lines are the result of a fit with simulated distribu- 
tions (dashed lines) corresponding to the population of dif- 
ferent vibrational levels T*2(v) and HD(u), respectively. The 
resulting population probabilities are shown in the inserts; for 
histogram bins spanning more than one unit in v, the relative 
populations were forced to be equal within the bin. 

Ref. [10]), the distance d= ((x 1 -x 2 ) 2 + (yi-y2) 2 ) 1/2 be- 
tween the two fragments in the detector plane is plotted 
rather than the weighted projected distance D defined by 
Eq. ©. The projected distance distributions P(d) carry 
information about the kinetic energy release as well as 
about the angular distribution of the fragments with re- 
spect to the ion (electron) beam direction [12]. Since 
at Ed = eV the angular correlations are expected to 
be isotropic for symmetry reasons, P(d) can be readily 
analyzed to yield the kinetic energy release distribution. 
In both channels only the vibrational states built 
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on the electronic ground state of the fragments 
are energetically accessible, i.e., the available final 
states are D 2 (X 1 E g (w))+H(ls) with u = to 20 and 
HD(X 1 E 9 (w))+D(ls) with v = to 17, respectively. The 
relative populations of these states, which are derived 
from a fit of the projected distance distributions P(d) 
using simulated distributions for individual vibrational 
states, are shown in the inserts. In the simulation we 
assume the rotational temperature of the parent ion to 
be 100 K (see also further below) and the DR rate coeffi- 
cients to be constant for all rotational angular momenta. 
For these fits some vibrational levels were grouped to- 
gether and an equal population was assumed within each 
group to reduce the uncertainties caused by the anticorre- 
lation of contributions between energetically close levels; 
they are nevertheless still dominating the accuracy of the 
deduced populations. 

The derived population distributions of the vibrational 
levels for the molecular fragment are similar for both 
channels and similar to the distributions measured by 
Strasser et al. HU in the DR of H^ and D^. They all 
show a bell shape-like behavior, close to what is expected 
from a simple phase space argument (To| . 



D. The three-body channel H+D+D 

Also in studies of DR fragmentation channels which re- 
sult in a total breakup of the molecule into its atomic con- 
stituents, the EMU system has some advantages as com- 
pared to standard imaging techniques based on micro- 
channel plates and an optical readout, even though the 
latter usually exhibits a considerably better position res- 
olution (e.g., < 100 /im [l2| as compared to ~ 750 fim of 
the EMU detector). Besides the unambiguous fragment 
identification, the advantage is in particular the high data 
rate of up to 2000 Hz that can be handled by the EMU 
system (to be compared to, e.g., ~ 30 Hz of the TSR 
optical imaging setup fl2T|). 

In the DR of D 2 H + at Ed = eV the three atomic 
products can only be formed in their electronic ground 
states. The total kinetic energy release Ek in the 
D+D+H channel is thus uniquely determined by the 
ground state energy release Qo of D 2 H + with respect to 
the asymptotic free neutral atomic products plus the re- 
maining excitation energy of the D 2 H + molecule, which 
we describe by a (rotational) temperature. As discussed 
in Sec. |Hj the weighted projected distance D 2 , which 
is not affected by any of the ambiguities in the frag- 
mentation geometry determination, is proportional to 
the transverse energy release. The distribution of the 
weighted projected distances D 2 thus provides informa- 
tion on the total kinetic energy release. 

The measured distribution, derived from data collected 
after 30 s of electron cooling, is plotted in Fig. [TT] to- 
gether with the simulated distribution, which was calcu- 
lated assuming that the three fragments are isotropically 
distributed in the available phase space and also assum- 
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FIG. 11: Normalized distribution of weighted transverse dis- 
tances D 2 observed in the DR of D2H + at Ed = eV for 
the fragment channel H+D+D after 30 s of phase-space cool- 
ing. The solid line is the result of a simulation assuming an 
isotropic breakup geometry and an initial rotational excita- 
tion of the D2H + ions corresponding to a rotational temper- 
ature of 100 K. 



ing Qq = 4.67 eV and a rotational temperature of 100 K. 
While the edge of the distribution at large D 2 , which is 
most sensitive to the total kinetic energy release, is very 
well described by the simulation assuming a subthermal 
rotational temperature of 100 K (see also Ref. [||), there 
are small but statistically relevant deviations on the left 
shoulder of the distribution, which are likely caused by 
anisotropies in the momentum distribution between the 
three fragments. 

Although we only have access to the transverse kinetic 
energies of the fragments, it has been shown before @, 
[llj | that information about the momentum distribution 
between the three fragments can be gained employing a 
slightly modified concept of Dalitz coordinates and Dalitz 
plots. The coordinates r\\ and ?y 2 originally introduced by 
Dalitz [f| are linear combinations of the kinetic energies 
Ei of the three fragments in the centcr-of-mass frame, 
taking into account energy and momentum conservation. 
By plotting the number of events as a function of rji and 
?; 2 , the momentum correlations between the fragments 
can be visualized. In particular, for a purely statistical, 
phase-space dominated breakup the Dalitz plot is evenly 
populated. The modified concept [f| [TJ [id], [lj| consists 
in defining projected Dalitz coordinates Qi and Q 2 by 
replacing Ei by the corresponding transverse energy E±^ 
and by plotting the number of events as a function of Q\ 
and Q2 instead of 771 and r\i- The measured distribution 
is then divided by a simulated distribution assuming a 
phase-space dominated breakup in order to account for 
possible detection efficiency variations and to regain a 
uniformly populated Dalitz plot in case no momentum 
correlations between the fragments exist. We will refer 
to these plots as Dalitz ratio plots. 

In the case of the DR of D 2 H + where two of the three 
fragments are identical, the projected Dalitz coordinates 
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FIG. 12: Dalitz ratio plot of the three-body fragmentation 
channel following the DR of D2H + at Ed = after 30 s of 
electron cooling. The triangles depict projected dissociation 
geometries for a sample of points in the (Q\,Q2) plane; the 
white dot shows the position of the H, the black points the 
positions of the Ds. 



Qi and Q2 are conveniently denned as [U 




with 



E 



(8) 
(9) 

(10) 



The allowed values of Qi and Q2 are 
2 



and E±=J2 E± }i 
confined by Q\ + Q 2 2 <l/9 

In filling the Dalitz ratio plot two ambiguities have 
to be considered. The first is connected with the in- 
distinguishability of the two Ds, which leads to two en- 
tries: at (Qi,Q2) and (— Qi,Qa)- The second one is 
specific to the EMU system and is caused by the ambi- 
guity [see Figs, rjja) and Hfb) and Eq. (fl"U|) ] in deducing 
the breakup geometry from the hit pattern. This leads 
to two additional entries at (Q' lf Q2) and (— Q' x , Q2) and 
to a smearing of structures along the Qi coordinate. The 
resulting Dalitz ratio plot, divided by the distribution for 
an isotropic breakup of D2H + , is shown in Fig. Q21 It is 
quite obvious that the momentum distribution is far from 
being isotropic and that linear decay geometries are en- 
hanced by up to 40%, where one of the Ds and the H are 
emitted back-to-back with the second slow D remaining 
in the center. 

The present result agrees with the result obtained in 
Ref. using standard 2D-imaging, but the almost two 
order of magnitude higher event rates that can be han- 



dled by the EMU system are resulting in improved statis- 
tics and thus in a considerably clearer picture. The statis- 
tics which can now be attained in these measurements 
should finally be sufficient to employ the Monte Carlo 
image restoration technique discussed in Ref. with 
strongly reduced artificial noise patches. This will al- 
low for a more quantitative analysis of the momentum 
distributions. 



VI. SUMMARY 

The potential of the EMU imaging system for study- 
ing the dissociative recombination of polyatomic molec- 
ular ions with electrons in merged beam experiments is 
clearly borne out by the results obtained for the DR of 
D2H + . The main advantage of the new set up, which 
is based on a large area, energy and position sensitive 
Si detector with multi-hit capabilities, is the possibility 
to determine the individual masses of the fragments on 
an event-by-event basis. This allows one to efficiently 
distinguish DR events leading to only neutral fragments 
from DR and background events involving charged frag- 
ments, to uniquely identify these fragmentation channels, 
and to determine their branching ratios as a function of 
the relative ion-electron energy. While the position res- 
olution, which is limited by the width of the read-out 
strips to ~ 750 fim, cannot compete with the resolution 
of < 100 yum of optical imaging systems, it is neverthe- 
less sufficient to allow for detailed studies of total kinetic 
energy releases and transverse breakup geometries of the 
different fragmentation channels by 2D imaging; the lack 
of position resolution is at least partly compensated by 
the ~ 100 times higher event rate that can be handled 
by the EMU system, resulting in data of high statistical 
quality. Moreover, in DR studies of polyatomic molec- 
ular ions involving heavier fragments, also the smaller 
minimum distance between two fragments that can be re- 
solved by the EMU detector as compared to MCP/CCD- 
based 2D-imaging systems can be an additional advan- 
tage. 

Comparing the results obtained in the present investi- 
gation of the DR of D2H + with results of earlier measure- 
ments, where available, generally good agreement is ob- 
served. In particular, the branching ratios of 78.0(0.4) %, 
13.5(0.3) % and 8.6(0.2) % measured at E d = eV for 
the D+D+H, HD+D, and D2+H channels, respectively, 
confirm the results recently obtained at CRYRING Q 
using the transmission-grid method. The branching ra- 
tios as a function of collision energy were measured for 
the first time. While the general trend of the ratio of 
the combined two-body to the three-body fragmentation 
channel looks rather similar to what has been observed 
before for HT i and H 2 D+ @, the relative branching 
ratios between the D2+H and HD+D channels display 
an unexpectedly large isotope effect. Whereas the ratio 
2B(D 2 +H)/B(HD+D) is expected to be 1 in the absence 
of any isotope effect, this ratio is found to be enhanced 
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by about 25% at most energies investigated, and more- 
over reaches a so far unexplained large value of 3.7(0.5) 
at E d ~ 5 eV. 
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